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Abstract 

We analyze the single transverse-spin asymmetry (SSAs) for inclusive hadron production in the pp collision 
at RHIC based on the twist-3 mechanism in QCD. As an origin of SSAs, we take into account of all kinds of 
pole contributions associated with the twist-3 quark-gluon correlation functions in the polarized proton. By 
the inclusion of the soft-fermion-pole contribution in addition to the soft-gluon-pole contribution, we find SSAs 
observed at RHIC are better described for all kinds of mesons. Py-dependence of the asymmetry and the 
comparison with the FNAL data are also presented. 
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1 Introduction 

Since the first observation of large single-spin asymmetry (SSA) in p^p — > irX by the FNAL-E704 
Collaboration [I] , SSA has been one of the hot topics in high energy spin physics. Large SSA has also 
been reported in semi-inclusive deep inelastic scattering (SIDIS), ep^ —> ehX (h — ir, i^)[2|3j. In 
the pp collision, BNL-RHIC reported a similar large SSA at even higher energies [4] [5] [6j [7] . These 
large SSAs can not be explained within the conventional framework based on the parton model and 
perturbative QCD [8] which had been successful in describing various hard inclusive processes. They 
reflect novel quark-gluon substructure of hadrons beyond the parton model picture. 

By now QCD mechanisms leading to such large SSA have been understood in some detail, and 
they arc often classified into two categories depending on the kinematic regions they cover. One is 
the so-called "T-odd" distribution and fragmentation functions in the framework of the transverse 
momentum dependent (TMD) factorization. This mechanism describes the SSAs in the region of the 
small transverse momentum of the final hadron as a result of a correlation between the spin and the 
intrinsic transverse momentum of partons [9j [TO] [TTJ [TSJ H3 HH [15] . Although this approach provides 
a physical interpretation of SSA, proof for the TMD factorization has been limited to the almost- 
back-to-back jets production in e + e~ annihilation |16j. Drell-Yan |17j and SIDIS 18 . In addition, the 
universality property of the TMD distribution function becomes unclear, in particular, for hadron- 
hadron collisions [TO] . From the phenomenological analyses of SIDIS data given by CERN and DESY, 
a parametrization of the responsible TMD functions has been obtained [2D] . 

Another mechanism describes the SSA as a twist-3 effect in the framework of the collinear factor- 
ization [HI [221I231 EH ESJESl E71 EH EllEO]. This mechanism is designed for the description of SSA 
in the region where the transverse momentum of the final hadron can be regarded as a hard scale in 
the process. In this mechanism twist-3 quark-gluon correlation functions replace the parton distribu- 
tion and/or fragmentation functions, and gives the information on the multiparton correlations in the 
hadron. Although the above two mechanisms describe SSA in different kinematic regions, recent studies 
have shown that for a certain component contributing to SSA they give the equivalent description for 
SSA in the intermediate region of the transverse momentum for which both frameworks are valid, and 
these constitute a unique and unified QCD description for SSA [51] . 

In this paper we study SSA for the hadron (n, K) productions in the pp collision at the RHIC 
energies (yflS = 200 and 62.4 GeV), 

pt(p,Sx)+p(j/)->HPh) + x, (i) 

in the framework of the collinear factorization, extending the previous works [25] [24] [26] [30] . Here S± 
represents the spin vector for the transversely polarized nucleon, and p, p' , Ph are the 4- momenta of 
each particle. The SSA for the above process is defined as An = (cr^ — <j^)/(cA + a^) where <j^> are the 
cross section for the process ([T]) corresponding to S± (—S±). In this framework, three kinds of twist-3 
cross section contribute to the single-spin-dependent cross section Act = a 1 * — a^: 

Act = G (3) (xi,x 2 )«) f{x')®D(z)®a A 
+ h{x)® f(x')®E {z) (z 1 ,z 2 )®a B 
+ h(x)®EW(x' l ,x' 2 )®D(z)®a c , (2) 

where (8> denotes convolution with respect to the light-cone momentum fraction, and /, h, D are, 
respectively, twist-2 unpolarized distribution, transversity distribution and the fragmentation function. 
Two variable functions with the upper index (3) are the twist-3 multiparton correlation functions: 
G^(xi,X2), E^(x' 1: x' 2 ) and E^(zi,z 2 ) are, respectively, the distribution function in the polarized 
nucleon, distribution function in the unpolarized nucleon and the twist-3 fragmentation function for 
the final hadron. a a , <Jb , &c are the perturbatively calculable partonic hard cross sections. It has 
been shown in [24] that the third term in ([2]) can be neglected compared to the first term because of 
the smallness of the partonic cross section; be <C <JA- The second term of © involves the twist-3 
fragmentation function, for which one needs different formalism to calculate the cross section [32] . In 
this circumstance, phenomenological analysis of the data has been performed including only the first 
term of ([2]) 23, 26 . According to the twist-3 formalism for SSA coming from the twist-3 distribution, the 
cross section occurs from a pole part of an internal propagator in the hard part (See e.g. [22] [25] for the 



detail of the formalism) . These poles fix one of the momentum fractions in the twist-3 distribution, and 
those poles are classified into the soft-gluon-pole (SGP) and the soft-fermion-pole (SFP) corresponding, 
respectively, to the zero momentum of the gluonline and the quarkline in the twist-3 distribution. A 
characteristic feature of the SGP contribution is that some of them accompany a "derivative" of the SGP 
function, which enhances the cross section compared with the usual "nonderivative" terms. Accordingly 
it has been considered to be a main source of SSA compared with the SFP contribution, and the 
previous analyses of An data focussed on the former [231 0S3 ■ Even though the SFP contribution does 
not appear with the derivative, there is no clue that the corresponding nonperturbative function itself 
is small compared with the SGP function and therefore its importance depends on the magnitude of 
the partonic hard cross section. In a recent paper [30 , one of the authors has shown that the partonic 
hard cross section for the SFP contribution is much larger than those for the SGP contribution, and the 
total SFP contribution could be a large source of SSA if the SFP function has a comparable magnitude 
as the SGP function. 

The purpose of this paper is to present a numerical analysis of An for the light-hadron production at 
RHIC, p^p — > hX (h = 7r, A"), including all the effects from the twist-3 quark-gluon correlation functions 
contributing to the first term of ((2|). In principle, independent twist-3 quark-gluon distribution functions 
should be determined by performing a global analysis of various SSA data for SIDIS and the pp collisions. 
However, SSA data obtained at HERMES and COMPASS are in the low Pt region, and therefore is not 
suited for the analysis by the twist-3 formalism. We thus focus on the asymmetry in the pp collision in 
this analysis. Previous analysis of A n in the pp collision included only the SGP contribution [23l [26] . 
and could reproduce the An obtained at FNAL and RHIC reasonably well, albeit with a deviation 
for A^ and the Pt dependence J351 IS] ■ We will see in this paper that the SGP contribution itself is 
insufficient to understand the patterns shown in the RHIC An data for the n and K production and 
inclusion of the SFP contribution significantly improves the agreement with the data. This suggests 
that the SFP contribution could be an important candidate for the origin of SSA. 

The rest of the paper is organized as follows: In Sec 2, we summarize the necessary cross section 
formula as well as the twist-3 quark-gluon correlation functions. In Sec 3, we present our method of 
numerical analysis and the result. Section 4 is devoted to a brief summary. 

2 Cross-section formula 

There are two independent twist-3 quark-gluon correlation functions, Gf and Gf, contributing to the 
first term of ([2]). They are defined from the nucleon matrix element of the composite quark-gluon 
operator on the light-cone (See [351130] for the detail): 

Z7T J ZTT j 

= ^(the^^GFix^xz) + i^-iTsfiijSlGWxuXi) + ■■■■ (3) 

where n is a lightlike vector satisfying p ■ n = 1 , ipf is the quark field for flavor a with the spinor index 
i, F al3 is the gluon's field strength, e apnS± = e a ^Pp,j,n v Sxp (com = +1), and • • • denotes twist-4 or 
higher. The nucleon mass Mn is introduced to define Gf and Gf dimensionless. From parity and T 
invariance, G F {x\,X 2 ) and G F (xi,x 2 ) satisfy the relation 

G F {x 1 ,x 2 ) = G a F {x 2 ,x 1 ), G F (x 1 ,x 2 ) = -G a F (x 2 ,x 1 ). (4) 

The twist-3 correlation function for antiquark flavor can be defined similarly to (j3j , replacing the nonlocal 
operator in ((3|), ^pj(0)gF al3 (fin)nfj^f(Xn), by its charge conjugated operator. The "antiquark" twist-3 
correlation function obtained this way is related to the original "quark" correlation function as 

G F (xi,x 2 ) = G a F (-x 2 ,-xi), G F (xi,x 2 ) = -G a F (-x 2 ,-x l ). (5) 

The contribution to the single-spin-dependent cross section from the correlation function ([3]) has been 
obtained in the literature. For completeness we give the explicit formula here in our notation: The SGP 
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where a s = g 2 /(4ir) is the strong coupling constant, the Mandelstam variables for the process are defined 
as 

s = (p + p') 2 , 
T = ( P -P h ) 2 , 
U = (p'-P h ) 2 , 



and the lower limits of the integration are given by z min = —^—^-,x' min — s+l j/ z ■ fa( x ) and D a (z) 
are, respectively, distribution and fragmentation functions for quark (or antiquark) flavor a. G{x) and 
D g (z) are, respectively, distribution and fragmentation functions for gluon. Partonic hard cross section 
in ([5]) are expressed in terms of the Mandelstam variables in the parton level: 
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where N c — 3 is the number of colors for SU(3). In ©, G F does not contribute to the SGP contribution 
because G F (x, x) = due to the symmetry property ((4]), and G F (x, x) appears in the form of — '"'"' 
G a F (x,x) [261 [28]. The SFP contribution was obtained in [30] as 



dx 



d*Ao^P a 2 s M N n MSi 



pO 

h dP 3 



dz f 1 dx' f dx 1 . / -x'U/z 

o \ x — 



S 2 J Zmm z 3 J x , mm x' J x x'S + T/z V x'S + T/z 

J2(G a F (0,x)+G F (0,x)^ x 

~a,b,c 

{f b (x')(D c (z)a^ c + D- c (z)a F b ^) + h{x')(D c {z)o F al ^ c + D- c (z)a F ^ 5 )} 
+ J2 ( G f M + G a F (0,x)) (f b (x')D g (z)a^ g + f h (x')D g (z)a^ g ) 



a . 6 



+ Y, ( G fM + G a F (0,x)) G(x')(D c (z)a F g ^ c + D- c (z)a F g ^) 

a,c 

+ ^(g£(0,z) + G£(0,x))g(x')A,(^ 



T F 

ag^g 



where the partonic hard cross sections are given by 



^ F 

a ab^c 



(N 2 s + 2i)(s 2 + u 2 ) -{N 2 t + u - s)s 



N 2 t 3 i 



NHi 



8ab5 a 



Oab^c = °> 



~F 
"ab^c 



~F 
"ab-tc 



a ab^g 



~F 
"ab^-g 



(N 2 c u + 2t){s 2 +u 2 ) . (iV 2 u + 2s)(P+u 2 h (i\r 2 -l)n 2 . 

"»c "I , T „^„ r ^ O a b r — ; ~„ a bd ac , 



N 2 t 3 U 

-(N 2 t + 2s)(t 2 +u 2 ) 



<lib 



N 2 s 2 tu 
-N 2 s + i-u 



N 3 st 2 



N 3 u 2 



VabOac-) 



•l 



+ — tz^(N- :(s 3 + 3s 2 u - 2u 3 ) + s 6 - s 2 u)5 ab 



+ 



N 2 s 2 tu 
{N 2 s + 2t)(s 2 +u 2 ) 

N 2 i 3 u ' N 3 siu 2 

-(N 2 u + 2t)(s 2 + u 2 ) 
N 2 i 3 u 
1 / u 1\ 1 fs 2 + si+i 2 u\ ( N c (u 3 -i 3 )(i 2 + u 2 ) 
N 3 {~§i + ti) + N < 



dab 



ag^-c 



SU 2 t 2 / STM' 

N 2 (s 3 - ii 3 )(s 2 + u 2 ) su(s 2 + su- it 2 ) - N 2 (s 4 + s 3 ii + ,s 2 m 2 + .to 3 + u 4 ) 



ag—^c 



{N 2 - l)st 3 u 2 
(N 2 u + 2s)(P+u 2 ) 
N C (N 2 - l)s 2 iu ' 
S + 2t -N 2 s - (N 2 t + 2s) (t 2 + u 2 ) 

J Vac H~ 



7V 2 (^ 2 -l)si 2 M 2 



N 2 (N 2 - l)u 2 



N C (N 2 - \)s 2 tu 



(14) 

(15) 

(16) 
(17) 

(18) 

(19) 

(20) 

(21) 
(22) 



-N 2 



(N 2 - l)s 2 t 3 ii 2 

1 



N 2 (N 2 - l)st 2 u- 



(4s 6 + llsH + 19sH 2 + 22s 3 P + 19s 2 £ 4 + list 5 + U 6 ) 



■{-stu z + N 2 (s 4 + sH + 2s 2 P + st 6 + P)) 



(23) 



The summation ^ a b , ^ a b , ^ a in © and (fT4")l indicate that the sum of a is over all quark and 
antiquark flavors, and the sum of b and c is over all quark flavors for the quark flavor a and is over all 
antiquark flavors for the antiquark flavor a. As is shown in (U4[) . the SFP contribution appears in the 
form of G F (0, x) + G F (0,x), which allows us to parametrize G F (0,x) + G F (0,x) when we study the 
process p^p — > hX. 

For completeness we also list the twist- 2 unpolarized cross section for pp — > hX [33] which becomes 
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d 3 a unpo1 _ o| f 1 dz f 1 dx' C dx I ( ~x'U/z 



dP 3 h 



S 



x' J x x'S + T/z \ x'S + T/z 



a,b,c 

+f a (x)f b (x')(D c (z)a^ c + D 5 (z)a^ 5 )} 

a 

]T fa(x)G(x')(D a (z)a^ a + D g (z)cr u ag ^ g ) 

a 

£ G{x)f a (x'){D a (z)Sr u ga ^ a + D g {z)a u ga ^ g ) 

a 

]T G(x)G(x')(D g (z)a^ g + D a {z)a u gg ^ a 



(24) 



' ab— >c 



' ah— >c 



'ab-yc 



"ab-tc 



rr U - 

~U 
O ag^ra 

±.U 

ag^rg 

~U 
® ga^ra 

±.U 

°ga^g 
fV 

gg^a 
a U 

gg^g 



C F s 2 + u 2 

n c i 2 ' 



C F s 2 + g 

N r U 2 ' 



2C F s 2 

n 2 iii 



5 ac 5, 



be- 



C F P 


+ u 2 


N c 


~o ^a 
S z 


C F P 


+ u 2 x 


N c 


s z 


2C F 


I 2 + p 


N c 


lit 


c F 


u 2 + s 2 


N c 


us 


c F 


s 2 + P 



C F s + u 

~N~c t 2 
C F s 2 + P 
N c u 2 



Snr - 



'Jbc 



zu F w 

2C F P 



^ac^bci 



N 2 su 



QacObci 



2G 



u 2 + P 



u 2 + s 

p 



N c st 

.u 
ag^g 

.U 
ag^ta 

1 ii 2 +P 



1 u 2 + P 



2N C ut 2C F s 2 



9 / ut us st 



where C F = (N 2 



2 \~ s 2 p u 2 
1)/(2N C ) with N c — 3. With the above defined cross section, An for p^p 
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3 Numerical analysis 

Independent unknown functions appearing in An are the SGP functions 

G a F (x,x), (37) 

and the SFP functions 

G F (Q,x)+G a F (0,x), (38) 

for each quark and antiquark flavor (a = u, d, s, u, d, s). Until now, the An data obtained at RHIC 
are the only ones to which twist-3 mechanism in the collinear factorization can be directly applied. In 
fact, the next-to-leading-order (NLO) QCD can describe the unpolarized cross section for the pp — > irX 
process at RHIC energy, which constitutes the denominator of An- We thus use the An data by the 
STAR Collaboration at ^/S = 200 GeV gj [5] and those by the BRAHMS Collaboration at y/S = 62.4 
GeV [7]. All data for A N from the STAR Collaboration are in the range of PhT > 1 GeV, so we adopt 
all of the data presented in (4) [5] in the fitting. On the other hand, some of the BRAHMS data are 
in the low PhT region. But they reported data for -n^ and K , which are important for disentangling 
the flavor structure of the quark-gluon correlations in the nucleon. So we will adopt data for 7r and 
if in [7] with PhT > 0.7 GeV in the fitting for this initial study on the SFP contribution. The total 
number of data points turn out to be 61. Unlike the analysis in [26], we shall not use FNAL-E704 data 
in the fitting, since the NLO QCD can not reproduce the unpolarized cross section at the fixed target 
energy [31] . 

For the parametrization of the SGP and SFP functions, we assume the following form: 

G F (x,x) = N G x a °(l - xf°f a (x), (39) 

G a F {0,x) + G a F (0,x) = N^x a '(l-x) fi 'f a (x), (40) 

where f a {x) is the unpolarized quark distribution, and N^,a^,/3^, N^,a^,f3^ are the parameters to 
be determined by the x 2 fitting to the data. Of course, the twist-3 quark-gluon correlation functions 
Gp(xi,X2, ju) and G F (xi,X2,n) have their own logarithmic scale dependence (/i dependence) which 
differs from that of the twist-2 unpolarized parton densitie^ and therefore the above ansatz for the 
SGP and SFP functions is not correct in the rigorous sense. However, at the present stage the SSA 
data is not so accurate, as the precise scale dependence of the twist-3 correlation is required. We 
therefore assume that the SGP and SFP component of the twist-3 correlation functions obey the same 
scale dependence as the twist-2 distribution through the ansatz (|59")1 and (|40p for simplicity. For the 
unpolarized parton densities in the nucleon, we use the LO parametrization of GRV98 [36]. As for the 
fragmentation functions for pions and kaons, we use the parametrization in |37) (hereafter referred to as 
the DSS fragmentation function), which can describe a greater variety of data in wider kinematic ranges, 
in particular for kaons, compared with the older more common ones in the literature. When we calculate 
An we set /i = P^t for the STAR data and p = 1 GeV for the BRAHMS data in all distribution and 
fragmentation functions. In the % 2 -fitting analysis below, we have used MINUIT package program in 
the CERN library. 

Since our main purpose here is to investigate the role and importance of the SFP contribution in 
comparison to the SGP contribution, we will analyze the data in terms of the following three sets of the 
fitting functions: 

FIT 1: SGP + SFP contributions with 13 free parameters: 

N G N? N G = N G N G = N G 



1 For the complete LO evolution equation for Gf(ii,I2,c) an d Gf(xi,X2, fi), see |35| . 
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Table 1: Result for the fitting parameters for FIT 1 (SGP+SFP: 13 parameters), FIT 2 (SGP only: 
Seven parameters) and FIT 3 (SGP only: 12 parameters). 



u - = n u R u — R u R u - = R u 

u,d,8,u,d,s — " ' ^u,d — h'v ■> ru,d,8,s ~ ^s 

I F N? N F = N F N F = N? 

^ -a F , (3 F , - H -- = P F . (41) 



FIT 2: Only SGP contribution with seven free parameters: 



N G , JVf, N? = N°, N<* = N£, 

a°, -,- = a G « G -« G « G - <<' 



#f, d = #f, ft AatS = ft. (42) 



FIT 3: Only SGP contribution with 12 free parameters: 

N G N? N G N G N- 

„G „,G „,G _ „,G oG oG oG _ oG 



<> <#> < a , d - s - = «,: ft, K, ft^ s = ft. (43) 

Several comments are in order for the above constraints among the parameters: With the present 
limitation of the accumulated SSA data, we could not obtain a convergent set of parameters when we 
treated all N G ' F , a G ' F and (3 F ' G as independent in the fit with SGP+SFP contributions. We thus 
introduced the above constraint relation in FIT 1. By allowing reasonable different values for N G ' F , we 
expect to be able to assume a G ' F and (3 G ' F are relatively flavor-independent. This is particularly so for 
a G ' F , since An rises in the forward direction and is expected not to be very sensitive to a F ' G . On the 
other hand, SGP appears as derivative contribution and thus we assume different values for /3 G between 
the valence flavors (a — u, d) and the sea flavors (a — u, J, s, s). As for the relations among N F ' G for 
the "sea" flavors, we have tried all possible combinations among N a ' — ±JVr ' and N^' — ±AT g ' . 
It turned out that the above constraint relation in FIT 1 gives the best fit. 

We have also performed FIT 2 and FIT 3 with only the SGP contribution in order to see the role of 
the SFP contribution. FIT 2 assumes the same relation as the SGP contribution in FIT 1 but the SFP 
contribution is omitted, i.e. FIT 2 has less degrees of freedom compared with FIT 1 by the lack of the 
SFP contribution. In FIT 3 with the SGP contribution, we allowed more freedom in the parameters. 
This fitting contains the Fit II of [26] (hereafter, simply referred to as KQVY06) as a special case as 
well as the above FIT 2. With this maximally allowed freedom in the SGP contribution, the number 
of degrees of freedom in FIT 3 is almost the same as FIT 1 and thus one can see whether the SFP 
contribution can be replaced by an appropriate choice of the SGP contribution. 

The results of the fitting in the above three cases are shown in Table 1. The calculated An arc 
also shown in Fig. 1 together with the RHIC data. For comparison, we have also shown in the same 
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Figure 1: The results of the three fits for ^4^. The solid lines are for FIT 1, the long-dotted lines are 
for FIT 2, and the dashed lines are for FIT 3. The result by KQVY06 is also shown by the short-dotted 
lines. 
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each quark-flavor obtained in FIT 1 in comparison to the unpolarized quark distribution f a (x) [scaled 
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Figure 3: SGP functions at Q 2 = 2.4 GeV 2 obtained in three fits. The SGP function in KQVY06 is also 
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the sign of this relation is opposite. 



figure the result of KQVY06 [2B] obtained with the CTEQ5L distribution function [55J and the Kretzer's 
fragmentation function [39] . From Fig. 1, one sees that FIT 1 and FIT 3 reproduce A jv data for all 
mesons, while FIT 2 and KQVY06 fails for K~ . We note that our FIT 3 has converged into a result 
that is different from KQVY06. This is because we have used different sets of data from KQVY06. 
The deviation of the curves of KQVY06 from data for the charged mesons is also due to the difference 
of the used data sets. (Also note the kinematics of [7 differs from Fig. 10 of 26 .) This feature is 
also represented in the resulting x 2 , which are 61.47, 132.9 and 64.90 for FIT 1, FIT 2, and FIT 3, 
respectively. Although FIT 3 reproduces An equally well as FIT 1, the behavior of the obtained fitting 
functions are extreme and looks unphysical as is seen by the resulting parameters in the functions shown 
in Table 1. This suggests that the SFP contribution is also a necessary source for SSA and plays an 
important role which cannot be substituted by the SGP contribution. 

Figure 2 shows the obtained SGP and SFP functions in FIT 1 for each quark flavor in comparison 
with the unpolarized parton density (scaled by factor 1/10). One sees that the SGP and SFP functions 
are comparable in magnitude, but the SFP functions are more shifted in smaller x regions. Among 
SGP functions, those for the "valence" flavors (u and d quarks) spread more into a larger x region than 
those for sea flavors. This spreading into a larger x region actually plays an important role for the 
rising An in the forward region. We observe that the flavor dependence of the SFP functions also has a 
similar tendency. Owing to the behaviors of each function, the largest contribution to A N and A N are, 
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respectively, from u and d quark SGP functions. As for the kaons, the SGP contribution from strange 
quark turned out to be the largest with small corrections from other flavors and the SFP contributions. 
This is largely due to the large strangeness component in the DSS fragmentation function. From the 
definition of the SGP and SFP functions, the latter might be interpreted as a more "exotic" component 
of the correlations, while the former may be closer to a quark parton density accompanied by the soft 
gluon. This may lead to the SFP function to have nonzero values only in the smaller x region. 

To compare the functions in FIT 1 with other cases, we have shown in Fig. 3 the SGP functions 
obtained in each fit. For comparison we have also shown the SGP functions of KQVY06. As is seen from 
Fig. 3, the SGP functions for u and d quarks in FIT 3 are huge compared with those for other quark- 
flavors, and the former is sharply peaked, which are the result of the extreme values of the parameters in 
FIT 3 shown in Table 1. This unphysical behavior of the SGP function can be taken as a manifestation 
of the fact that the SGP function alone can not reproduce all of the observed An , and other sources of 
SSA are necessary. From Fig. 3, the SGP functions in FIT 1 turn out to be smaller compared with those 
in FIT 2 and FIT 3, since the SFP functions can give rise to the additional ^4jvs. Moderate behaviors of 
the SGP and SFP functions given in FIT 1 as shown in Figs. 2 and 3 indicate that the SFP contribution 
is a natural source of An together with the SGP contribution. 

To see the relative magnitude of the SFP contributions in FIT 1, we have shown in Fig. 4 the 
decomposition of An into the SGP and SFP contributions. One sees from the figure that in many cases 
the SGP brings a larger contribution as a whole, while the SFP works as a small but non- negligible 
correction. For the kaons, we observed that the SFP contribution from u and d quarks are quite large 
but they tend to cancel each other. As one can see from © and (fT4"|) the SFP hard cross sections are 
larger than the SGP hard cross section typically by factor N% [3D]. However, the SGP functions spread 
more in the larger x region as shown in Fig. 2 and also contribute as a derivative. Therefore it brings 
a large contribution to An in the forward region. 

Figure 5 shows the P hT dependence of A N for v^ = 200 GeV together with the STAR data. All fits 
give similar curves for each xf and approximately reproduce experimental data. This is quite natural 
because our three fits used correct PhT at each xf of the An data shown in Fig. 1. If we had included 
E704 data in the fit without taking into account of the PhT dependence, deviation from the RHIC data 
may have been observed in the P^T-depndence as was the case in the analysis by KQVY06. 

Using the SFP and SGP functions in the fits, we have calculated A N at the FNAL energy y/S = 20 
GeV for p^p —>■ ttX and p<p — > irX. As mentioned before, we did not include the An data of FNAL-E704 
(y/S = 20 GeV) in the fitting. We do this comparison as a semiquantitative test of our fit result in 
terms of the twist-3 formalism at a lower energy. With our convention for the twist-3 distributions for 
the "antiquark" flavor in (JS|), those in the antiproton are related to those in the proton as 

Gp(Xi, X2) | anti-proton = Gp (X\ , X2) \ proton, G p{X\ , £2) | anti-proton == Gp{Xl, X2)\ proton- (44) 

The result was shown in Fig. 6 together with the E704 data pQ. From Fig. 6, one sees that FIT 1 and 
FIT 2 give qualitatively similar curves as the data, while FIT 3 shows the deviation from the data at 
smal\-XF for tt°. As pointed out in [33], the NLO QCD in the collinear factorization can not reproduce 
the unpolarized cross section at the fixed target energy and thus our twist-3 formula for An may not be 
directly used for the FNAL data. With this in mind, KQVY06 included the E704 data in their fitting by 
multiplying the ad hoc factor -/Ve704 = 1/2 to the theory calculation. If we had known this iVE704 more 
precisely, we could have included the E704 data in our analysis, which is beyond the scope of the present 
study. From Fig. 6, one sees the E704 data is smaller than the calculated An, which is consistent with 
the analysis in |26j . 

4 Summary 

In this paper, we have presented a numerical analysis of the RHIC An data for p^p — > hX (h = n, K) 
in terms of the twist-3 mechanism in the collinear factorization. As a source of the asymmetry we 
have included all the contribution from the twist-3 quark-gluon correlation functions in the transversely 
polarized nucleon, which consist of the SGP and SFP contributions. We have shown that the combination 
of the SGP and SFP contributions can reproduce the An data at RHIC with moderate choice of the 
functions, while it is difficult to reproduce all of them if one keeps only the SGP contribution. In 
particular, An for K was properly reproduced by the inclusion of the SFP contribution. Although the 
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relevance of the SFP contribution is numerically supported, we do not have any dynamical information 
on the shape and magnitude of the SGP and SFP functions at this stage. □ In addition, our present 
analysis did not include the second term of ([2]), which could be another important source of SSA. To 
clarify these points, one needs more variety of SSA data together with the nonperturbative study on 
the SGP and SFP functions. 
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There is, however, a model calculation I40| showing the SFP function is small. 
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